Localization dynamics of fluids in random confinement 
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The dynamics of two-dimensional fluids confined within a random matrix of obstacles is investigated using 
both colloidal experiments and molecular dynamics simulations. By varying fluid and matrix area fractions in 
the experiment, we find delocalized tracer particle dynamics at small matrix area fractions and localized motion 
of the tracers at high matrix area fractions. In the delocalized region, the dynamics is diffusive at long times, 
while in the localized regime anomalous diffusion at intermediate times and trapping in finite pockets of the 
matrix is observed. These observations are found to agree with the simulation of an ideal gas confined in a 
weakly correlated matrix. 



Introduction. An understanding of the dynamics in disor- 
dered heterogeneous media is of great interest for a wide range 
of fields like materials science, geophysics or biology HKD. 
The mass transport in such media is associated with a strong 
separation of time scales, i.e. there are typically mobile parti- 
cles that move through a "matrix of immobile" particles. Ex- 
perimental realizations of such systems are binary mixtures of 
colloids with disparate sizes (|4] |5], other examples are ion- 
conducting glasses @{TOl or biological systems like cells f3). 

To model such heterogeneous systems one often considers 
the matrix particles as being fixed in space. One of the sim- 
plest models of this type is the disordered Lorentz gas lfmiT2l 
where a single tracer particle moves in a matrix of fixed, over- 
lapping hard spheres that are randomly distributed in space. 
This model exhibits a delocalization-to-localization transition 
which is associated with a critical percolation point [13 1 of 
the void space accessible for the tracer particles. Whereas at 
low matrix densities the void space forms a percolating net- 
work leading to diffusive motion of the tracer particle, above 
the critical percolation density the void space becomes discon- 
nected and the tracer particle is trapped in finite pockets of the 
matrix. At the critical point of the localization transition, the 
long-time diffusion of the tracer particle is anomalous lfT4l - 
[161, as reflected by a sublinear increase of the mean-squared 
displacement (msd), 5r^(f) ^ f*^ with x w 2/3 in two dimen- 
sions (2D) and jc w 0.32 in 3D lH]. 

It is an open question as to whether the localization scenario 
of the Lorentz model is also relevant in more realistic and 
complex heterogenous media like ion-conductors or porous 
materials. Simulations have reported anomalous diffusion and 
localization dynamics in a range of more realistic heteroge- 
neous systems such as binary mixtures with a frozen-in matrix 
component lfT7l - E2l or with a disparate size ratio of the species 
[23 1 as well as realistic models of ion-conducting alkali sili- 
cate glasses (see, e.g., [81). The latter studies are corroborated 
by calculations in the framework of the mode coupling theory 
of the glass transition ||24| . This theory provides a qualita- 



tive description of the Lorentz gas 11251 and predicts localiza- 
tion transitions as in the Lorentz gas for disparate-sized binary 
mixtures [26| and quenched-annealed systems Il27l430l . 

While models like the Lorentz gas are very instructive for 
the understanding of the dynamics in heterogeneous media, 
experimental realizations of systems with fixed, randomly- 
placed matrix particles are very rare. Most of the experimen- 
tal studies of particle diffusion in disordered media have been 
reported for porous glasses fT\\ where systematically control- 
ling and varying the pore size distribution and connectivity is 
very difficult. In this Letter, we present a 2D colloidal model 
experiment consisting of small tracer particles in a random 
confining matrix of fixed large particles, that shows delocal- 
ized tracer particle motion at low matrix area fractions and 
highly localized tracer dynamics at high matrix area fractions. 
The experimental observation of localized dynamics is sup- 
ported by molecular dynamics (MD) simulations where we 
consider a 2D binary mixture of soft spheres with the large 
species being fixed and the small species forming an ideal 
fluid of non-interacting particles. 




1.5 



CO 



1.0 



0.5 





Line 
— — Line 1 

— Line 2 

— Simulation 



Bt H 



10 



20 



Figure 1: a) State diagram for the effective area fractions of the 
fluid (^f) versus the matrix particles ("I'm)- b) Schematic of a two- 
dimensional cell with the 4.95 |J.m matrix particles acting as spacers, 
c) Structure factors for the matrix-matrix contributions for lines 0, 1 
and 2 and the simulation. 
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Experiment. A binary colloidal mixture of 3.9 |i.m and 
4.95 |J.m diameter super-paramagnetic polystyrene spheres 
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(microParticles) dispersed in water is confined between two 
glass slides to make a two-dimensional (2D) sample cell as il- 
lustrated in Fig. [TJi. The large particles act as spacers between 
the glass slides and hence form the fixed matrix in which the 
small particles - the fluid - are free to move |,32il . Video mi- 
croscopy is used to image an area of 430 [a.m x 340 |J.m for 
typically up to two hours. Standard particle tracking proce- 
dures are used to find the matrix and fluid particle coordinates 
||33l as a function of time from which we computed the static 
and dynamics quantities of interest. 

We control the effective area fractions of both the ma- 
trix (^m) and fluid (^) particles in-situ by applying a per- 
pendicular external magnetic field B which leads to a re- 
pulsive pair potential given by UpMir) — — ^ where r 
is the distance between two particles, jUq the permeability 
of free space and Xf.m the magnetic susceptibility of the 
fluid or matrix particles. To obtain (pu and (pp we calcu- 
late the fluid and matrix particle effective hard sphere diam- 
eters (7f,m using a Barker-Henderson approach according to 
(ypM = <M + .Co^ (^1 _e-/3£/F.M('-)^ dr, where a^^ are the 
hard sphere diameters in the absence of a magnetic field and 
j3 = l/ksT ll34l l35l . This approach enables us to very effi- 
ciently probe the dynamics of the fluid particles at different 
effective matrix and fluid area fractions without changing the 
matrix configuration. Note that the resultant effective matrix 
to fluid particle size ratio remains fairly constant at 0.787 as 
is also evident from the linear paths in the (0f, ^m) state dia- 
gram (Fig.[T^). Our experimental state points are represented 
by three lines, denoted line 0, 1 and 2, in the (0m, state 
diagram as shown in Fig. [T^. The lowest point in each line 
was obtained with no magnetic field. In our analysis, we will 
focus on the lowest and highest state points of line 1 and 2, 
which are labeled as 'LlPl' and 'L1P6' for line 1 and 'L2P1' 
and 'L2P4' for line 2, respectively. 

To characterize the structural correlations within 
the matrix we computed its static structure factor 
= ^<Jlf^lJ|,.e-'^ (^'-oOy Here A^m is the num- 
ber of matrix particles, (.) represents an ensemble average 
and r, and r, are the positions of matrix particles / and /', 
respectively. The structure factors for the matrices cor- 
responding to line 0, 1 and 2 only show weak fluid-like 
structural correlation (Fig. [T];), similar to the uncorrected 
Lorentz model. Note that the structure factors for the matrices 
do not change along each line, which is a direct consequence 
of our approach to increasing the effective matrix and fluid 
area fractions by increasing the magnetic field. Importantly, it 
also directly confirms the fixed nature of the matrix particles. 

We analyzed the fluid particle dynamics of all state 
points by computing the corresponding MSD according to 

5r2(f) = ^<^£^j[ri(f)-ri(0)]2^, with A^f the number of 

fluid particles and r, the position of fluid particle /. The re- 
sulting MSDs along lines 1 and 2 are shown in Fig. [2^. Note 
that the shift in the short-time diffusion between the two lines 
is due to the different number densities It36i 37 1. At LlPl, 



characterized by low matrix and low fluid area fractions, the 
tracer dynamics is diffusive at long times. Upon increasing the 
effective area fractions 0f and 0m along line 1 , the dynamics 
show a noticeable slowing down and strongly sub-diffusive 
behavior at intermediate times. However, at long times the 
MSDs of all the state points along line 1 exhibit diffusive be- 
haviour, corresponding to delocalized motion. Consistently, at 
all state points along line 0, i.e. at lower matrix area fractions, 
the tracer particles show delocalized motion as well (data not 
shown). Remarkably, the fluid particle dynamics exhibit com- 
pletely different behavior along line 2, i.e. at much higher 
matrix area fractions. At L2P1, the MSD exhibits a plateau at 
long times, which upon increasing the effective area fractions 
(j)F and along line 2, decreases to smaller values. This be- 
haviour clearly indicates the localization of the tracer particles 
in this region of the state diagram and shows that the length 
scale associated with the localization becomes smaller as the 
matrix area fraction is increased. 

Simulation. To gain more insight into the mechanism of 
the experimentally observed localized tracer dynamics, we 
performed detailed MD simulations of a quenched-annealed 
system: a binary system of purely repulsive soft spheres 
with WCA interaction potential, where the matrix particles are 
fixed. In particular, we simulate the tracer dynamics in the 
limit of (j)p 0, while we systematically increase the matrix 
area fraction, which enables us to efficiently sample the re- 
gions corresponding to both delocalized and localized tracer 
dynamics. 

In the MD simulations, the matrix particle diameters are 
sampled equidistantly, Om G [0.85, 1.15], to avoid crystalliza- 
tion. The fluid particles are non-interacting to avoid the ef- 
fect of fluid-fluid interaction on the dynamics. The interaction 
potential is Uapir) = 4e„|j[(cT„^/r)l2 - ((7„^/r)'^] + for 

r < '"cut = 2^^^<yap and zero otherwise, with a,/3 G (F,M). 
The particle diameters are given by Gap = {Oa + Op)/!, 
the energy scales by Emm = 1, Emf = 0.1, and Eff = and 
the masses are set to my^ = mp = 1. Continuity of energy 
and forces at rcut is ensured with the smoothened potential 
Uapir) X (r - r,ut)V[h^ + {r - r,ut)^] with h = 0.005(Jmm. 
The velocity- Verlet algorithm is used to integrate Newton's 
equations of motion with time step St = 7.2 x lO^'^fo, where 

f0 = [mM(TM/(eMM)]^''^ = l- 

For each realization of the system, between 1000 and 
4000 matrix particles are equilibrated at number density p = 
0.278(7j^^ at the temperature A:Br = 1 by randomly selecting 
the particle velocities from the Maxwell distribution every 100 
steps, for 10^ time steps. Then, the matrix particles are fixed 
and their positions are uniformly rescaled to p = 0.625c7jy[^. 
The computed quantities are averaged over 100 independent 
matrix configurations. To simulate the tracer dynamics we in- 
sert up to 1000 non-interacting fluid particles into each matrix 
configuration. We vary the tracer particle diameter, which is 
equivalent to varying the matrix area fraction without chang- 
ing the structure of the matrix, analogous to the experiments. 
At each diameter the system is equilibrated for 10^ time steps, 
after which we start our production runs in the microcanon- 
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Figure 2: a) MSDs for the fluid particles in the experiment along lines 1 and 2. b) Mean square displacements for the fluid particles in the 
simulation for fluid particle diameters Op = 0.2, 0.3, 0.4, 0.45, 0.5, 0.6, 0.7, and 0.9. The dotted line roughly marks the limit above which 
finite size effects become important (with simulation box size L). 
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Figure 3: a) Self-part of the van Hove correlation function for the fluid particles along line 1 at state points 1 and 2 after 250s and 3000s. b) 
As a), but for line 2. Self-part of the van Hove correlation function for the fluid particles in the simulation for c) a delocalized state, d) an 
intermediate case and c) a localized state. 



ical ensemble. Varying the system size L allows to control 
finite size effects, which only start to play a role after the MSD 
exceeds as indicated in Fig.|2j3. 

Comparison between simulation and experiment. We com- 
puted the structure factor for the matrix particles in the simu- 
lation, which is shown in Fig.[T];). Consistent with the exper- 
iments, it shows weak fluid-like structural correlations. 

In Fig.|2|3 we present the mean square displacements for the 
simulations for increasing CTp - which is equivalent to increas- 
ing 0M- In the long-time limit, the MSD changes from diffu- 
sive at small diameters - i.e. low ^ - to localized at large di- 
ameters - i.e. high 0m - with an extended subdiffusive regime 
at intermediate times. This clearly indicates that the system 
undergoes a transition from delocalized to localized dynam- 
ics. While at short times the motion of the tracer particles is 
ballistic (5r^ °^ t^) and thus different from the experiment, the 
curves qualitatively match the experimental data in Fig. at 
intermediate and long times. The anomalous diffusion is most 
pronounced for 0.4 < Op < 0.5, where we find a power law. 



8r^ ^ t^, with xPil/l) the exponent predicted for the Lorentz 



model 1 161 . The rounding of the localization transition in our 
simulations is a direct consequence of the soft sphere inter- 
actions |38|. The simulated MSD curves closely match the 
experimental data in Fig.|2^: the experimental MSDs for line 1 
closely resemble the simulated MSDs in the delocalized region 
(Op < 0.4), while the MSDs for line 2 are consistent with the 
simulated MSDs in the localized region (Cp > 0.5). 

We also computed the self-part of the van Hove correlation 
function ||39l Gs(r,f) for the experiments and the simulations 
as it gives the full spatial information for the tracer particle dy- 
namics at a given time. This function gives the distribution of 
displacements Ar,- = — n(0) | of a tagged particle / at time 



3 1, we show the ex- 



f,Gs(r,f) = 4(irii5(r-Ar,)).InFig. 

peiimental van Hove coiTelation functions Tor LlPl and L1P6 
at times t = 250 s and t = 3000 s. The lower time (250 s, solid 
lines) roughly corresponds to the time where sub-diffusion is 
strongest, while the larger time (3000 s, dashed lines), marks 
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the end of the experimental runs. The delocahzed nature of 
the fluid particles at low area fractions is reflected by the broad 
distribution at LlPl and its time dependence: the width of this 
distribution increases by about a factor of 3 from t = 250 s to 
t — 3000 s. At point L1P6, the shape of Gs (r, t) is very similar 
compared to LlPl but the distributions are less broad at both 
times. The van Hove correlation functions for L2P1 and L2P4 
are presented in Fig. [3]3. While L2P1 is characterized by a 
far lower fluid area fraction compared to L1P6, the Gs(r,f) of 
these two points are fairly similar The localization of the sys- 
tem at L2P1 can only be inferred from the time -dependence, 
i.e. from the MSD. The Gs{r,t) for L2P4 indicates a strong 
localization as is inferred from the very narrow distribution 
and the absence of any shift when changing from t — 250 s to 
f = 3000 s. 

In the simulation, we display the van Hove function at three 
times to compare to the experiment, where the values of the 
MSD roughly agree with those of the experiment: t « 3fo, 
t « 39fo and an additional time that is order of magnitude 
larger, t « 550fo- The shape of Gs{r,t) of the simulation at 
cTf = 0.2 (Fig.[3]:), matches the Gs(r,f) of LlPl (Fig.[3}i) well. 
As c7f increases, Gs{r,t) in the simulation undergoes the same 
qualitative development as in the experiment. On the local- 
ized side at Op = 0.9 (Fig. [3]l), a very narrow distribution is 
observed that shows virtually no change of Gs(r, f) with time, 
except from a small broadening which qualitatively matches 
with our observations at L2P1 and L2P4. 

The general agreement between experiment and simulation 
shows that fluid-fluid interactions are not important in the ex- 
periment. This is also indicated by the fact that there is no sign 
of particle hopping in the experimental Gs{r,t) which would 
be reflected in peaks or shoulders around the nearest-neighbor 
distance between tracer particles (note that also the quantity 
2;rrGs(r,f) does not show such a feature) ISl [TSl l38l . 

Conclusion. A combination of an experimental colloidal 
model system and MD simulation has been used to study the 
dynamics of two-dimensional fluids confined in a random ma- 
trix of obstacles. In the experiments, mean squared displace- 
ments and van Hove correlation functions show signatures of 
delocalized tracer dynamics at low matrix area fractions and 
localized motion at high matrix area fraction. In particular, 
we observe long-time diffusion at small matrix area fractions, 
while anomalous diffusion at intermediate times and trapping 
in finite pockets of the matrix is present at high matrix area 
fractions. The simulation of an ideal gas confined in a weakly 
correlated matrix shows qualitative agreement with the exper- 
iment. Around the percolation point of the localization tran- 
sition, the mean squared displacements, as obtained from the 
simulation, exhibit power-law-like anomalous diffusion over 
multiple time decades compatible with the predicted exponent 
of the two-dimensional Lorentz model. 

Our experimental model system provides an ideal system 
for the investigation of transport in strongly heterogeneous 
media. The interactions between the tracer particles are tun- 
able, the number density of the tracer particles can be varied 
and it is possible to introduce external fields by using optical 



tweezers or, in the case of charged tracer particles, an elec- 
tric field. Thus, it is possible to study the change of transport 
properties from Lorentz gas-like localization scenario to the 
"usual" glass transition in porous confinement if one increases 
the density of tracer particles, and, in doing so, one can probe 
the nonlinear response to external fields. 
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